Phenomenology of a-axis and 6-axis charge dynamics from microwave spectroscopy of 

highly ordered YBa 2 Cu 3 6 50 and YBa 2 Cu 3 6 993 
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Extensive measurements of the microwave conductivity of highly pure and oxygen-ordered 
YBa2Cu306+i; single crystals have been performed as a means of exploring the intrinsic charge 
dynamics of a d-wave superconductor. Broadband and fixed- frequency microwave apparatus to- 
gether provide a very clear picture of the electrodynamics of the superconducting condensate and 
its thermally excited nodal quasiparticles. The measurements reveal the existence of very long-lived 
excitations deep in the superconducting state, as evidenced by sharp cusp-like conductivity spectra 
with widths that fall well within our experimental bandwidth. We present a phenomenological model 
of the microwave conductivity that captures the physics of energy-dependent quasiparticle dynamics 
in a d-wave superconductor which, in turn, allows us to examine the scattering rate and oscillator 
strength of the thermally excited quasiparticles as functions of temperature. Our results are in close 
agreement with the Ferrell-Glover-Tinkham sum rule, giving confidence in both our experiments 
and the phenomenological model. Separate experiments for currents along the a and b directions 
of detwinned crystals allow us to isolate the role of the CuO chain layers in YBa2Cu306+ M , and a 
model is presented that incorporates both one-dimensional conduction from the chain electrons and 
two-dimensional transport associated with the Cu02 plane layers. 



I. INTRODUCTION 

Microwave conductivity measurements have provided a 
powerful method for investigating the low energy charge 
dynamics of high temperature cuprate superconductors. 
Early on, the observation of a linear temperature de- 
pendence of the magnetic penetration depth X(T) in 
YBa 2 Cu306.95 provided some of the first evidence of d- 
wave pairing in these materials^ Over a decade later, 
the d-wave superconducting phase remains one of the 
cornerstones of our understanding of the cuprate prob- 
lem. Surprisingly, the physics of the superconducting 
state seems to be well described as a BCS instability of 
a two-dimensional Fermi liquid. However, it is essential 
that this conjecture be thoroughly tested, and one of the 
best ways to do this is through detailed measurements of 
charge transport in high quality single crystals. In this 
work we report on microwave spectroscopy of two gold- 
standard cuprate materials — Ortho-I YBa2Cu3 06.993 
and Ortho-II YBa2Cu3O6.50 — and develop a simple phe- 
nomenology that describes their low-energy charge dy- 
namics. 

Understanding the role of impurities in the cuprate su- 
perconductors has taken many years of great effort. The 
materials have complicated chemistry — containing at 
least four elements — and have much greater sensitiv- 
ity to disorder than conventional superconductors due 
to hnite-angular-momentum Cooper pairing. Early mea- 
surements on the temperature dependence of the mag- 



netic penetration depth of YBa2Cu3 06+j/ showed devia- 
tions from the linear temperature dependence expected 
for a clean d-wave superconductor. It was later suggested 
that these effects could be the result of pair-breaking 
disorder a conjecture that was confirmed by showing 
that the controlled addition of Zn impurities to high qual- 
ity YBa 2 Cu306+ y crystals caused a crossover in AA(T) 
from T to T 2 behaviour4i£i& Ongoing efforts to uncover 
the intrinsic physics of the cuprates depend on continued 
improvements in sample quality. The samples used in 
this work are the most well-ordered cuprates available, 
and are of such exceedingly high purity and crystallinity 
that even the off-plane dopant oxygen atoms are highly- 
ordered 

Some of the earliest measurements of the low en- 
ergy electromagnetic response of the cuprates were made 
by Nuss et al, who measured the conductivity of 
YBa2Cu307 thin films using terahertz spectroscopy over 
the frequency range 0.5 THz to 2.5 THz,— The observa- 
tion of a peak in u\ (T) below the critical temperature T c 
was attributed to a strong decrease in inelastic scatter- 
ing, rather than the conductivity coherence peak that oc- 
curs in s-wave BCS superconductors. Shortly afterward, 
Bonn et al. measured a similar peak in o\ (T) in the mi- 
crowave range ^Li2i They proposed that this was the re- 
sult of a competition between the rapid increase in quasi- 
particle lifetime r and the decaying quasiparticle density 
n n (T) as the temperature was reduced. This analysis was 
carried out within the context of a Drude model char- 
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acterized by an energy-independent transport relaxation 
rate r _1 . Other transport measurements including ther- 
mal conductivity*^ and higher frequency electromagnetic 
measurements^ have also confirmed the rapid collapse in 
scattering on entering the superconducting state. 

In a more recent publication, our group presented 
a series of microwave measurements that confirmed 
and extended these earlier findings^ This was made 
possible by two experimental developments: the ad- 
vent of higher purity YBa2Cu306+i, crystals grown 
in BaZrC>3 crucibles and the development of a suffi- 
cient number of fixed-frequency cavity perturbation sys- 
tems to map out a coarse microwave conductivity spec- 
trum. Five superconducting resonators, operating be- 
tween 1 GHz and 75 GHz, were used to examine the 
surface impedance of the same single crystal sample of 
fully-doped YBa2Cu3 06.993 over the temperature range 
2 K to 100 K (the conductivity inferred from these data is 
shown in Fig.[|)]of the present article). At each tempera- 
ture the conductivity spectrum was fit to a simple Drude 
model, which captured the features of the data well, thus 
supporting the assumptions used in the earlier work. The 
temperature dependent quasiparticle scattering rate r _1 
derived from these fits was seen to decrease rapidly below 
T c and saturate at a constant value of 6xl0 10 s~ x below 
20 K, corresponding to a spectral width of approximately 
9 GHz. Assuming a Fermi velocity Dp = 2 x 10 7 cm/s, 16 
the residual scattering rate implies a quasiparticle mean 
free path of 4 /zm, a very large distance for a complicated 
quaternary oxide. 

However, the observation of a temperature indepen- 
dent transport relaxation rate is at odds with stan- 
dard models for the electrical conductivity of a d-wave 
superconducto»i£ and resolving this elastic scattering 
problem is a major goal of the work presented in this arti- 
cle. The strong temperature dependence of the transport 
relaxation rate at higher temperatures, where the scat- 
tering processes are inelastic, has motivated models in 
which the opening of a gap in the spectrum of electronic 
excitations in the system leads to a collapse of the inelas- 
tic transport relaxation rateii^iiS Finally, the same set of 
fixed frequency resonators was used to probe the con- 
ductivity for currents along the crystal 6-axis of Ortho- 
I YBa2Cu3 6.993- In these measurements, currents are 
parallel to the CuO chains and reveal an anisotropy in 
charge conduction that can be attributed to the sum of 
a narrow Cu02 plane quasiparticle spectrum and a very 
broad spectral feature arising from the one-dimensional 
CuO chains^ 

Greater insight into the low temperature charge dy- 
namics of YBa2Cu306+j, was gained following the devel- 
opment of a novel low temperature microwave spectrom- 
eter that could measure the conductivity spectrum as 
a continuous function of frequencyi2i*22i Spectroscopy of 
Ortho-II YBa2Cu306.5o revealed that the nodal quasi- 
particle spectrum did in fact contain most of the fea- 
tures expected for weak-limit impurity scattering in a 
clean d-wave superconductor, namely a sharp cusp-like 



spectrum and a linear-in-temperature spectral width. In 
contrast, similar spectroscopy of Ortho-I YBa2Cu3 06.993 
material displayed much more Drude-like spectra, ex- 
cept at the very lowest temperatures. We return to 
examining these features in the present article, present- 
ing a complete body of microwave data on Ortho-I and 
Ortho-II YBa2Cu306+ y , and show how these data can 
be described by a simple phenomenlogy that captures 
the physics of energy-dependent charge dynamics. 

The paper is organized as follows. Section [n] begins 
with a brief discussion of the two main experimental 
techniques, fixed-frequency microwave cavity perturba- 
tion and broadband bolometric surface resistance spec- 
troscopy. In Sec. IIVI we use the broadband bolometric 
surface resistance data to extract detailed conductivity 
spectra and present a simple phcnomcnological model 
of the conductivity. With this in hand, in Sec. we 
process data from the five fixed- frequency experiments, 
which provide a means of capturing the conductivity 
spectrum at higher temperatures than can be reached 
with the bolometric experiment. Finally, we examine the 
two key parameters that emerge from fitting to the ex- 
perimental data: the width of the conductivity spectrum, 
which measures the thermally-averaged transport relax- 
ation rate (Sec. IV Aj l and the total integrated spectral 
weight or 'quasiparticle oscillator strength' (Sec. IVBfl . 
Throughout this article we use recently revised absolute 
values of the T — ► penetration depths obtained using 
a newly developed zero-field ESR absorption method^ 
The revised penetration depths are considerably shorter 
than existing values in the literature and result in sub- 
stantial quantitative changes to the previously published 
quasiparticle conductivity and supcrfluid density. How- 
ever, we emphasize that neither the qualitative form of 
the conductivity data nor the main conclusions drawn 
from it are affected. 



II. EXPERIMENTAL TECHNIQUES 

In this paper we focus on the microwave conductivity, 
er = or — i<72, in the long wavelength limit and in the 
linear response regime. In experiments that probe the 
conductivity, the quantity that is inferred from measure- 
ments is the microwave surface impedance Z s = R s +iX s . 
When the applied microwave fields are screened on a 
length scale that is large compared to the characteris- 
tic size of the wavefunctions of both the Cooper pairs 
and the single-particle excitations (the coherence length 
£ and mean free path £, respectively) it is possible to de- 
scribe the charge response by local electrodynamics, in 
which case: 

Zs = R s+ iX s = ^. (1) 

In very clean cuprate samples the mean free path will 
exceed the penetration depth and the electrodynamics 
will not necessarily remain in the local limit. However, 
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it is possible to take advantage of the two dimensional 
electronic structure of the cuprates, which gives rise to 
a short coherence length and mean free path along the 
crystal c direction, to ensure that the charge response is 
locali22i2^ In the work reported here we have exclusively 
used such geometries. 

As pairs condense in a superconductor, the reactive 
response of the superfluid quickly overwhelms the dissi- 
pative response of the quasiparticle excitations, and over 
most of the temperature range cr 2 ^S> <J\ . The local elec- 
trodynamic relation then simplifies to the following ap- 
proximate forms: 



fl.(w,T) 
X,(u,T) 



i^ w 2 A 3 (T)CTi(W)T)) 

(j, u>\(T). 



(2) 



From these expressions it is clear that a measurement 
of the surface resistance R s contains information about 
the real part of the conductivity oi, while a measure- 
ment of the surface reactance X s is a direct probe of the 
penetration depth A. As with nearly all forms of spec- 
troscopy, technical limitations require one to work with 
incomplete spectroscopic information and subsequently 
to rely on causality to extract the quantities of interest. 
Our present work is no exception. However, in this ar- 
ticle we present a wide range of measurements that not 
only allow us to extract the frequency-dependent con- 
ductivity, but also provide a set of checks and balances 
through the oscillator-strength sum rule that let us verify 
our results in detail. 

In these measurements, a well-defined electromagnetic 
geometry must be used for three reasons: to separate the 
individual components of the conductivity tensor, to en- 
sure that demagnetization effects are well controlled, and 
to enforce local electrodynamics. A particularly clean ap- 
proach that has been widely used is to place the sample 
in a microwave enclosure near a position of high symme- 
try, where there is an electric node and the microwave 
magnetic field is quasi-homogeneous. 

Single crystal samples of many of the cuprate super- 
conductors grow naturally as platelets having a broad a-b 
plane crystal face and thin c-axis dimension. In this case, 
demagnetization effects are minimized if the broad face 
of the sample is aligned parallel to the microwave mag- 
netic field. In response to the microwave fields, screening 
currents flow near the surface of the sample, across the 
broad a-b face, closing along the c direction to complete a 
loop. All measurements presented in this article employ 
thin crystals with negligible c-axis contribution, and were 
carried out in this low-demagnetizing factor orientation. 

An experimental feature common to each of the seven 
sets of apparatus employed in this work is the use of a 
sapphire hot finger m& This arrangement allows a sample 
to be held at an elevated temperature in the microwave 
fields, while the resonator remains at base temperature. 
For the field configuration we utilize, the sapphire plate 
is transparent to microwaves and makes a negligible con- 
tribution to the dissipation. The sample is affixed to 



the sapphire plate (with a very small amount of sili- 
cone grease) that is then inserted into the microwave 
field through a hole in the metal end-wall that is small 
enough to cut off microwave propagation at the operating 
frequency. The thermometer and heater are located out- 
side the cut-off hole, in intimate thermal contact with the 
sapphire isothermal stage, which is in turn weakly linked 
to the 1.2 K pumped helium bath. 

Resolving the microwave surface impedance of small 
single-crystal cuprate samples over a wide temperature 
range has been achieved by the use of cavity-perturbation 
techniques iSLSS In these experiments the sample under 
test is brought into the microwave fields of a high quality- 
factor (high-Q) resonant structure such as a supercon- 
ducting cavity or a loop-gap resonator. In cavity per- 
turbation, shifts in the resonant frequency / and quality 
factor Q are related to the surface impedance viai^ 
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= iVAZ s = -T(AX S - W s ). (3) 



The quantity 6(l/Q) = l/Q load (T) - l/Q unload , where 
Q load (T) is the quality factor of the resonator when it 
contains the sample at temperature T, and Q unload is the 
quality factor of the empty resonator. Sf — f load (T) — 
/ load (T = 1.2 K) is the shift in frequency of the loaded 
resonator on warming the sample from base temperature 
to temperature T. Similarly, AX S = X S (T) - X S (T = 
1.2 K). R s is the absolute surface resistance, and the 
geometric factor T depends on the dimensions of both 
the cavity and the sample, as well as the geometry of 
the fields that are being perturbed, and is determined 
empirically. 

The direct relationship in Eq. [31 assumes that the sam- 
ple is much thicker than the length scale of the electro- 
magnetic screening — either the superconducting pene- 
tration depth A or the normal-state skin depth. In situ- 
ations when the screening currents penetrate to a depth 
approaching the sample thickness, Eq. |3| must be modi- 
fied. For our Yl^CuaOg+j/ samples, which have a thick- 
ness w 10 fim, this becomes important only in a very nar- 
row temperature region just below T c , and in the normal 
state where the skin depth is large. To account for this 
effect, we use a model that solves the electromagnetic 
field equations in our one-dimensional geometry with the 
microwave field applied parallel to the broad face of the 
sampled The result is a more complete version of Eq. [3| 
that explicitly involves the sample thickness t: 



1 - 



taah(/c£/2) 

Kt/2 



(4) 



where V s is the sample volume, V c is the effective cavity 
volume, /3 is a constant set by the type of cavity and 
its field configuration, and k = ^/kj/io<r is the complex 
propagation constant for the electromagnetic field in the 
sample. For this expression, Sf is modified and is now 
the shift in frequency on inserting the sample into the 
empty resonator, Sf = / load (T) - / unload . The physical 
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quantities of interest, including R s , X s , A, o\ and (T2 are 
determined from measurements of Sf and 5{\/Q) and the 
relationships a = k 2 /(i^iquj) and Z s — i^ioco/n. 

While it is possible to calibrate R s measurements to 
very good accuracy, cavity-perturbation techniques do 
not easily yield an absolute value for A(T — > 0) since a 
measurement of the absolute penetration depth amounts 
to a comparison of the total volume of the sample to the 
volume of the sample minus the small field-penetrated 
region at the surface of the sample. It is impossible to 
measure the dimensions of a single crystal sample pre- 
cisely enough to achieve an absolute measure of the in- 
plane A this way. Recent Gd ESR measurements have 
provided a breakthrough in absolute penetration depth 
measurements by using the characteristic microwave ab- 
sorption of Gd ions lightly doped through the crystal as a 
magnetic tracer to measure the field-penetrated volume 
at the surface of the sample directly. 23 

The calibration procedure for the superconducting 
resonators relies on measurements of well-characterized 
PbSn samples of similar dimensions as the YBa2Cu30g+ a 
crystal, as well as runs without any sample at all. These 
measurements allow for the small corrections necessary 
to account for background losses and the finite volume 
of the crystal. Each resonator experiment must be cali- 
brated independently. As a result, there is up to a 10% 
systematic uncertainty in the overall scale factor at each 
frequency. Stochastic noise is typically much less than 
5% for measurements well below T c . 

In cavity perturbation experiments, the resonator is 
generally restricted to operate at a single fixed frequency, 
requiring the use of many separate sets of experimen- 
tal apparatus to capture a frequency spectrum. Further- 
more, the cavity perturbation method requires that the 
dissipation of the unknown sample be comparable to the 
dissipation of the cavity itself in order for the absorption 
to be measured with high precision. This is a strong re- 
quirement on the cavity when the sample is a high quality 
superconductor in the T — ► limit. The measurement of 
the residual absorption in superconductors is challenging 
at any frequency: in the case of infra-red spectroscopy the 
problem becomes that of measuring values of reflectance 
that are very close to unity. In microwave cavity pertur- 
bation measurements there is a similar problem of dis- 
cerning the small difference in cavity Q with the sample 
in and out. Both techniques have problems with calibra- 
tion and discerning small differences between two large 
numbers. 

These sensitivity problems are circumvented in a direct 
absorption measurement. We have recently developed 
an apparatus capable of obtaining continuous-frequency 
measurements of the absolute microwave surface resis- 
tance in low-loss single crystal samples, which has been 
described in detail elsewhere»2i The instrument employs 
a bolometric method of detection in which the sample 
of interest is exposed to a microwave frequency magnetic 
field H T { whose frequency can be varied, and the corre- 
sponding temperature rise is measured. This method re- 



lies on the fact that the power absorption in a microwave 
magnetic field is directly proportional to the surface re- 
sistance R s : 

P abs = Rs I H 2 { dS, (5) 
Js 

where H r { is the root-mean-square magnitude of the mag- 
netic field at the surface S. Measurements of the tem- 
perature rise of the sample as a function of frequency 
directly give R s (u>). To enhance rejection of spurious 
temperature variations, the rf power is amplitude mod- 
ulated at low frequency and the resulting temperature 
oscillations of the sample are detected synchronously. To 
provide a frequency-tunable microwave magnetic field, 
the sample is placed in a custom-made rectangular coax- 
ial transmission line that supports a TEM mode in which 
the H fields lie in the transverse plane and form closed 
loops around the centre conductor. The line is termi- 
nated by shorting the center conductor and outer con- 
ductor with a flat, metallic end wall, thus enforcing an 
electric field node at the end of the waveguide. The sam- 
ple is positioned with its flat face parallel and very close 
to the endwall such that it experiences spatially uniform 
fields over its dimensions. A second metallic reference 
sample of known absorption is placed in an electrically 
equivalent position that serves to calibrate the absolute 
field strength, which varies considerably across the fre- 
quency range due to standing waves in the transmission 
line. This technique has provided a very sensitive means 
of measuring the broadband surface resistance spectrum 
in superconducting crystals over the range 0.5 GHz to 
21 GHz, permitting the first exploration of the detailed 
low temperature conductivity line shape in a number 
of compoundsi 22 i 31 i 32 The conductivity spectrum a(ui) 
can be extracted very simply from R s (ui) using Eq. 
Our full analysis is slightly more detailed, allowing for 
the small frequency dependence of the penetration depth 
that comes from thermally excited quasiparticles. 

A very convincing verification of the results from multi- 
ple sets of experimental apparatus is provided by a com- 
parison of the broadband R s (uj,T) data with measure- 
ments of the same sample in five different superconduct- 
ing resonators. As seen in Fig.^ the agreement between 
the experiments is excellent. In each of the six experi- 
ments the absolute value of R s is determined to better 
than 10%, with the principle uncertainty coming from 
the calibration procedure. 

III. SAMPLES 

The bilayer cuprate compound YBa2Cu306+j, has the 
distinct advantage that it has a stoichiometric cation 
composition with a very low level of cross-substitution, 
less than one part in 10 4 . Tunable hole doping of the 
Cu02 planes is achieved by varying the oxygen content 
of CuO chain layers, which act as a reservoir for oxygen 
dopants. However, the chains add the complication of 
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FIG. 1: Comparison of surface resistance measurements made 
on the same sample of YBa2Cu3Cv50 using the broadband 
experiment (solid symbols) with those from five microwave 
resonators (open symbols). The agreement between methods 
is excellent. The data are plotted as R s (uj)/lo 2 to remove the 
frequency dependence associated with the superfluid screen- 
ing. 



an extra one-dimensional conduction channel that acts 
in parallel with the two-dimensional channel associated 
with the CuC>2 planes. Varying the oxygen content of 
YBa 2 Cu 3 06+y usually occurs at the expense of intro- 
ducing off-plane disorder. Fortunately there are sev- 
eral highly-ordered phases for the chain oxygen atoms 
that are stable at room temperature, and we exploit 
two of these here. The first is the fully-oxygenated 
YBa2Cu30g.993, with Ortho-I order in which every CuO 
chain is nearly completely filled and T c = 89 K. Opti- 
mal doping for Yl^CusOg+j, occurs at y = 6.93, corre- 
sponding to T c = 94.0 K. Increasing the oxygen content 
beyond 6.93 reduces T c . Since adding oxygen to the CuO 
chains is a diffusion-limited process, filling a millimetre- 
sized crystal beyond y = 6.993 is impractical. For this 
publication, we have combined measurements on two dif- 
ferent Ortho-I crystals to provide a complete data set: 
sample #1 has T c =89 K and sample #2 has T c =91 K. 
The two samples provided very similar results but with 
some interesting differences that we attribute to the small 
difference in oxygen contained in the CuO chains. 

The other ordered phase we have studied is Ortho-II, in 
which the CuO chains alternate between full and empty. 
Ortho-II order is seen to exist over a wide range of doping, 
from y = 6.35 to y = 6.60, however the most complete 
chain filling has been found to occur at y = The 
particular sample used for the present microwave study 
is very close to this, having y = 6.50 and T c = 56 K. 
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FIG. 2: The surface resistance Rs(T) for Ortho-II ordered 
YBa2Cu3Oe.50 measured with currents in the a-direction (top 
panel) and 6-direction (bottom panel). 



Figure [2 presents fixed-frequency microwave ab- 
sorption measurements for a high quality crystal of 
YBa2Cu306.5o made in a geometry where the screening 
currents are along either the a or the b crystal direction. 
The surface resistance is seen to drop rapidly, by up to 
four orders of magnitude, on cooling below T c , and dis- 
plays a strong frequency dependence in the microwave 
range. Much of the frequency dependence of R s is due to 
the factor of w 2 in Eq.|2lthat results from screening by the 
superfluid. The remaining frequency dependence is due 
to the quasiparticle conductivity o\ itself. The raw sur- 
face impedance data for fully doped YBa2Cu3 06.993 have 
already been published and are not reproduced herei^2£ 
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IV. PHENOMENOLOGY 

For a superconductor, the conductivity can be ex- 
pressed as the sum of two conduction mechanisms acting 
in parallel: one due to charge carriers in the supercon- 
ducting condensate (denoted by a s f) and another due to 
thermally excited quasiparticles (denoted by cr qp ): 



10 



<j(uj,T) = a si (uj,T) + a qp (uj,T). 



(6) 



The superfluid conductivity can be parameterized using 
the London phenomenology^ 



n s (T)e 2 



7t5(lo) 



(7) 



where the temperature dependent superfluid spectral 
weight is represented by n s (T)e 2 /m* . This is an experi- 
mentally accessible quantity that is related to the London 
penetration depth in the following manner: 



n s (T)e 2 



1 



rrr 



MoA 2 L (T) 



(8) 



The goal of this work is to explore the phenomenology 
of the superconducting state without assuming a micro- 
scopic model for the conductivity of thermally excited 
d-wave quasiparticles. We extract cr qp (a;,T) from our 
measurements of the real and imaginary components of 
Z s (uj,T) using Eq. ^ to obtain both the spectral line- 
shape and its temperature dependence. This process 
would be straightforward if we could measure both R s 
and X s at each frequency. However, measurements of 
the frequency dependence of X s are very difficult to carry 
out, for reasons that are closely related to the difficulty 
in measuring absolute penetration depth at a fixed fre- 
quency. To carry out the analysis of our broadband 
data, we therefore use measurements of X s at one low 
frequency only (1 GHz) and extend this to higher fre- 
quencies using a self-consistent analysis of the surface re- 
sistance data. We are helped greatly by the fact that the 
frequency dependence of A is very slight, due to the dom- 
inant contribution from the superfluid. The small fre- 
quency dependence of A is due solely to the weak screen- 
ing contributions from the thermally excited quasipar- 
ticles. Note that we distinguish between the magnetic 
penetration depth A = l/^/Mo^o^ and the London pene- 
tration depth A l that appears in Eq. 1121 



A. a-axis Phenomenology 
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FIG. 3: Broadband absorption measurements of R 3 for 
Ortho-II ordered YBa2Cu3O6.sc) with screening currents par- 
allel to the crystal a-axis (top panel) and 6-axis (bottom 
panel). Plotting the data as R s /f 2 gives a quick means of 
estimating the conductivity spectrum from the raw data (see 
Eq. Here we emphasize that the frequency dependence 
does not follow l// 2 as it would for a Drude conductivity, but 
rather the power is closer to 1/f over most of the range. 



We begin by analyzing the microwave conductivity for 
currents parallel to the crystal a-axis. This is perpen- 
dicular to the CuO chains, so these should make a neg- 
ligible contribution to the a-axis conductivity. In addi- 
tion, we expect that the Cu02 bilayer should lead to two 
nearly identical contributions to o~ a (u), since the splitting 
between bonding and antibonding combinations of the 



planar wavefunctions is expected to be weaki^ We will 
therefore model the a-axis conductivity by a single con- 
ductivity spectrum due to one effective two-dimensional 
band. 

There is currently no accepted phenomenology to de- 
scribe the detailed quasiparticle conductivity of a Cu02 
plane. Several author siS^SiSi have assumed a Drude-like 
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form: 



(w,T) 



n„(T)e 2 



l/r(T) 



(9) 



where n n (T)e 2 /m* is the temperature dependent quasi- 
particle spectral weight (often referred to as 'normal 
fluid') and l/r(T) is the thermally- weighted and energy- 
averaged quasiparticle transport relaxation rate. While 
this particular model has proven useful for interpreting 
qualitative features of the electrodynamics, there is no 
a priori reason why it should model the quasiparticle 
conductivity correctly. In fact, broadband measurements 
demonstrate that the lineshape is distinctly non-Drudc- 
likei— *2£ A very clear and simple demonstration of this 
is obtained by plotting R s / f 2 at fixed temperature T, 
as shown in figure Fig.|31for YBa 2 Cu306.5o- Within the 
low temperature approximation of Eq. [21 Rs/ f 2 is pro- 
portional to <ti(w). The key point is that R s / f 2 does not 
follow a l// 2 power law at high frequencies, as would 
a Drude spectrum. The data are much closer to fol- 
lowing a 1/f power law. Consequently, we have pro- 
posed a more realistic phenomenological lineshape to de- 
scribe the a-axis quasiparticle conductivity spectrum in 
YBa 2 Cu 3 



6+y 



.22 



<Ti a (u,T) 



{ujA{T)}y( T ) 



(10) 



where o~o(T) represents the zero frequency limit of the 
quasiparticle conductivity, A(T) the inverse spectral 
width, and y(T) the anomalous power. Note that this 
is equivalent to a Drude model if we set y = 2. 

For the spectrum with y ^ 2, several conditions apply 
to the use of Eq. EE 39 The expression applies to posi- 
tive frequencies only, and y must be greater than one or 
the oscillator-strength sum rule will not be satisfied (i.e. 
o~x(u>) will not be integrable). A more subtle point is 
that the exponent y must always eventually reach a high 
frequency value of 2 since the energy-dependent relax- 
ation rate is cut off by the thermal windowing. However, 
Eq. ^| provides excellent fits to the spectra presented in 
this article. Furthermore, Eq. 1101 is very useful in model- 
ing the quasiparticle contribution to the imaginary part 
of the conductivity, because it can be extrapolated and 
integrated beyond the frequency range of our measure- 
ments. It will be shown below that success in preserv- 
ing the oscillator strength sum rule in the superconduct- 
ing state supports the use of this phenomenological form 
as a means of accounting for the normal fluid contribu- 
tion to screening. Since the electrical conductivity is a 
causal response function, the real part of the conductiv- 
ity determines the imaginary part via a Kramers-Kronig 
relation^ 



2w, 



a 2a (u;,T) = V dto'^ 

IT I ., t.\' 



,a la {uj',T) 



lj' 2 ~ UJ 2 



(11) 



Using this in conjunction with the measured superfluid 
response (via Eq. [7J yields a complete expression for the 
imaginary part of cr(w, T) in Eq. Q 



(T 2 (w,T) 



+ a 2a {uJ,T). 



(12) 



where V denotes the principal part of the integral. 



The strategy is to fit broadband R s (u>,T) data at a 
given T using the phenomenological model for the quasi- 
particle conductivity (Eqs. EH an d lll|) along with mea- 
surements of the magnetic penetration depth to parame- 
terize the superfluid conductivity (Eq. UJ . For the spec- 
trum at each temperature, we obtain best-fit parameters 
for <t (T), A(T) and y(T). Finally, we generate exper- 
imental estimates for the real part of the quasiparticle 
conductivity by solving Eq. \\\for ai a (u),T) and assum- 
ing <J2a{u,T) as given by the best-fit parameters. In 
many respects this procedure is akin to that used in 
infrared spectroscopy studies of superconductors, where 
one needs to extrapolate the measured reflectance to zero 
frequency in order to utilize Kramers-Kronig relations. 41 
The difference is that we need a model lineshape to ef- 
fectively extrapolate our measurements beyond our max- 
imum measurement frequency / max , but we do so ulti- 
mately for the same purpose. This technique requires 
sufficient curvature within the experimental bandwidth 
(i.e. A -1 < J max ) for a reliable extrapolation. This 
approach will be best for low T data, where A -1 will 
be governed by elastic scattering of d-wave quasiparti- 
cles from very dilute defects in our high purity samples. 
As we turn to modeling the data at higher T, inelas- 
tic scattering processes ought to push A -1 well beyond 
/ max . Fortunately, as the a\ a (uj,T) spectrum broadens, 
the quasiparticle contribution to a 2 (uJ,T) becomes negli- 
gible throughout our bandwidth. In this case, one can set 
a 2a (u,T) = and extract o~i a (u),T) from R s (u,T) and 
Al (T) without requiring any of the corrections discussed 
above. 



B. &-axis Phenomenology 

Figurc[3]shows the broadband 6-axis surface resistance 
data for YBa 2 Cu3C-6.5o plotted as R s /f 2 versus /. While 
the a-axis data seem to asymptotically follow l/f y at 
high frequencies, the o-axis data show a change in y 
beyond 10 GHz. In a previous publication we argued 
that YBa 2 Cu306+j/ should have a quasi-one-dimensional 
band whose low energy excitations in the superconduct- 
ing state are primarily chain-like, and will therefore only 
contribute to charge conduction along the &-axisi2S The 
chain scattering rate is expected to be much larger than 
that for the Cu0 2 planes, and should therefore give rise 
to a very broad contribution to a±b(uJ,T) that can be 
approximated as a temperature dependent constant over 
the experimental bandwidth in our analysis. The change 
in y observed above 10 GHz in the 6-axis spectra, but 
not in the a-axis, is due to this very broad component. 
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In this case the &-axis microwave conductivity can then 
be expressed as a sum of contributions from the two- 
dimensional bands (a term identical in form to Eq. 110(1 
plus a constant attributed to the quasi-one-dimcnsional 
band: 

^,T) = TT] ^ Ffy+ ^(^0 J T). (13) 

In general, the parameters erg, A and y for 6-axis currents 
will not be required to match those used to fit d-axis data. 
We extract aib(uJ,T) from R s data using measurements 
of Xl and the modeled <72&(w, T) as described in the pre- 
vious section (a constant in Eg. 1131 will give no contribu- 
tion to (Jib)- The procedure will have limitations similar 
to those for the d-axis case, but the addition of a fourth 
parameter a\ D renders the analysis less reliable as A -1 
approaches the experimental bandwidth at higher T. 

V. QUASIPARTICLE CONDUCTIVITY 

We show the real part of the quasiparticle conductivity 
spectrum for both in-plane principle directions at rep- 
resentative temperatures for both dopings in Figs. 
as obtained from our surface impedance measurements 
using the procedures described above. The figures are 
comprised of data from the broadband experiment (open 
symbols) and the microwave resonators (solid symbols). 
We include results from two different Ortho-I samples, 
denoted sample #1(T C =89 K) and sample #2 (T c =91 K), 
because we did not have a complete set of data (both 
crystal directions, all sets of apparatus) for a single sam- 
ple. As will be shown below, the two samples provided 
very similar results but with some interesting differences 
that we attribute to the small difference in oxygen con- 
tained in the CuO chains. In Figs. the smooth lines 
represent the best fit lineshapes obtained using Eqs. ITUI 
and 1131 to describe the d and 6-axis results, respectively. 
At each temperature our model lineshape appears to re- 
liably capture the high frequency behaviour, however our 
expressions do not accurately follow the data at the very 
lowest frequencies, particularly at T = 1.3 K. Note that 
our model for the two-dimensional conductivity term is 
necessarily concave downwards at low frequency while 
the lowest temperature bolometry data clearly suggest a 
cusp. Nonetheless, we argue that our model lineshape 
provides a reliable means of extrapolating our measure- 
ments to high frequency and provides a reasonable mea- 
sure of the width and spectral weight of the conductivity 
peaks. 

A key point made in the preceding section is that the 
quasiparticle conductivity spectra we have measured are 
distinctly more cusp-like than the Lorcntzian line shape 
of a Drude model. In Fig.|S]we present the best-fit values 
of y{T) that were used in Eq. ^5]to describe the conduc- 
tivity data extracted from the bolometry data. Error 
bars indicate the standard deviation obtained from fit- 
ting to R s (u,T) as a function of lo at each given tem- 




I i i i i i i i i i 

5 10 15 20 25 

Frequency (GHz) 



FIG. 4: Representative curves of the quasiparticle conduc- 
tivity ai(ui,T) for Ortho-II ordered YBa2Cu306.so measured 
with currents in the a-direction (top panel) and 6-direction 
(bottom panel). The fits shown use the phenomenological 
model of Eq. 1101 The low temperature cusp-like shape and 
high frequency power law (~ with y ~ 1.5) are fea- 

tures expected for a d-wave superconductor in the presence 
of weak-limit impurity scattering, and the spectrum broad- 
ens quickly as other scattering mechanisms are introduced at 
higher temperatures. 



perature. Two key features are important to observe: 
First, y{T) appears to be crudely a : b isotropic for each 
doping. Second, y(T) displays a relatively weak temper- 
ature dependence characterized by a slow increase with 
increasing temperature, and appears to saturate at ap- 
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FIG. 5: Representative curves of the quasiparticle conductiv- 
ity o"i (oi, T) for sample #1 of Ortho-I ordered YBa2Cu3 06.993 
measured with currents in the a-direction (top panel) and 
6-direction (bottom panel). The solid lines are fits to the 
data with Eq. 1101 and Eq. 1131 for the a and b-axis respec- 
tively. At the lowest temperature the cusp-like spectrum is 
evident, however a more Drude-like spectrum is recovered at 
higher temperatures, consistent with our earlier findings for 
this Ortho-I ordered YBa2Cu3 06.993- 



proximately 1.75 ±0.10 and 1.45 ±0.05 for the overdoped 
and underdoped samples, respectively. 

We now turn to the problem of extracting the real part 
of the quasiparticle conductivity from our cavity pertur- 
bation experiments. In this case we only have measure- 
ments at discrete frequencies, from independently cali- 
brated experiments. To help constrain the fits when us- 
ing such a sparse set of frequencies, we choose to fix y at 
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FIG. 6: The low temperature quasiparticle conductivity 
<ti(lu,T) for sample #2 of Ortho-I ordered YBa2Cu3 06.993 
measured using the bolometric method with currents in the 
a-direction (top panel) and ^-direction (bottom panel). The 
spectrum is very similar to that of sample #1, although some 
differences in the fit parameters will be discussed. 



the asymptotic high-temperature values of Fig. |H1 before 
proceeding to extract the real part of the microwave con- 
ductivity from R s {uj,T). We note that changing y ulti- 
mately leads to negligible changes in the extracted values 
of <7i (w, T) within our experimental bandwidth. How- 
ever, we will show below that the choice of y has implica- 
tions for the integrated quasiparticle spectral weight and 
that the values obtained are consistent with the oscilla- 
tor strength sum rule4£ Representative spectra from the 
resonator measurements with fits are included in Figs. Eh 
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FIG. 7: The exponent y of Eqs. 1101 and 1131 obtained from 
fitting to broadband measurements of the absorption spec- 
trum for YBa2Cu3 06.993 and YBa2Cu306.so samples. The 
exponent y describes the high frequency decay of the quasi- 
particle conductivity curves shown in Figs. The dashed 
lines indicate the asymptotic values of 1.45 (Ortho-II) and 
1.75 (sample #1, Ortho-I) subsequently used for fitting to the 
higher temperature spectra obtained with the fixed-frequency 
resonators. Sample #2 of the Ortho-I material shows larger 
values of y, corresponding to broader <ti(lj) curves, which we 
attribute to a slightly increased level of oxygen disorder in 
that sample. Despite this, the value of y is still inconsistent 
with the Drude value of y = 2. 



El 

We show <ti(uj, T) extracted from all of the cavity per- 
turbation measurements in Figs. [H] and Note in each 
plot there is a sharp peak in the curves at T c due to su- 
perconducting fluctuations ) 43 i 44 in addition to the large, 
broad peak which shifts to higher temperatures as the 
measurement frequency is increased. These results can 
be interpreted in the same qualitative manner as in previ- 
ous work. 15 ! 20 Neglecting the sharp peaks at T c , cri(w, T) 
increases smoothly on cooling from the normal state into 
the superconducting state. The large increase in a\ (u>, T) 
below T c can be attributed to a rapid decline in the 
quasiparticle transport relaxation rate below T c , com- 
bined with an increase in the condensate fraction—^ 
The frequency dependence of <7i(u>, T) becomes apparent 
at our measurement frequencies for temperatures below 
55 K for the Ortho-I sample and 24 K for the Ortho-II 
sample. This is in accord with the idea that the peak 
in the temperature dependence comes from the develop- 
ment of long quasiparticle transport lifetimes. 

The degree to which the extraction of ui{lo) from 
R s (iv) is affected by uncertainties in A(T — > 0) is an im- 
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FIG. 8: The temperature evolution of the quasiparticle con- 
ductivity <Ji(T) for sample #1 of overdoped YBa2Cu3 06.993 
measured with currents in the a-direction (top panel) and in- 
direction (bottom panel) . This data has been processed using 
the new absolute values of A(T —* 0) given in Table □ result- 
ing in a modification of the absolute magnitude of the curves 
from their previously published— form. 



Phase 


To(K) 


6+y 


Aa(nm) 


Ab(nm) 


A c (nm) 


Ortho-I 


89 


6.995 


103 ±8 


80 ±5 


635 ± 50 


Ortho-II 


56 


6.50 


202 ± 22 


140 ± 28 


7500 ± 480 



TABLE I: Experimental values of the anisotropic magnetic 
penetration depth measured by a zero field ESR technique^- 
The absolute value of A sets the absolute value of the conduc- 
tivity. 
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FIG. 10: The spectral width parameter 1/A(T) for both in- 
plane orientations and both oxygen dopings obtained from 
modeling the quasiparticle conductivity spectra shown in 
Figs. |1JJS| with Eqs. 1101 and 1131 The open symbols are from 
fits to the cavity perturbation data and the solid symbols are 
from fits to the bolometry data. The upper panel details the 
low temperature behaviour. 



FIG. 9: The temperature evolution of the quasiparticle con- 
ductivity <7i(T) for Ortho-II ordered YBa2Cu306.so measured 
with currents in the a-direction (top panel) and 6-direction 
(bottom panel). 



portant issue that must be addressed. From Eq. [3 it is 
clear that o\ « R s /X 3 . A 20% uncertainty in the abso- 
lute value of A therefore produces almost a factor of two 
error in the conductivity although it does not affect qual- 
itative features in <7i(w,T). In previous work, we have 
taken our A(T — > 0) values from the literature, usually 
from infrared reflectivity experiments^ which are able to 
resolve the a—b plane anisotropy. Very recently we have 
developed an accurate means of determining A/, (T — > 0) 
from broadband zero-field ESR absorption spectroscopy 



of Gd 3+ ions in Gd x Yi_ x Ba2Cu 3 6 +j / i^ The low tem- 
perature limiting values of the London penetration depth 
are quoted in Table [I] for samples of the same quality 
as those studied in this article, and they differ signifi- 
cantly from the values used previously>i£i2& Note that we 
quote our values of Xl(T — > 0) in Tabled to better than 
10% accuracy based on the redundant measurements and 
self-consistency checks that were built into our analysis 
protocol for that workiS We have reprocessed the pre- 
viously publishedi£*2S a and 6-axis conductivity data for 
YBa2Cu30g.993 and graphed the updated <j\{oj,T) data 
in Fig. 00 
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A. Spectral Width 

Having extracted the real part of the quasiparticle con- 
ductivity from all of our surface impedance measure- 
ments, we now turn to a detailed discussion of the best-fit 
parameters arising from the model in Eas. 1101 and 1131 for 
the d-axis and fr-axis data respectively. To begin, we show 
the spectral width parameter A _1 (T) for both in-plane 
orientations and both oxygen dopings in Fig. ^3 Repre- 
sentative error bars have been included at select temper- 
atures and indicate the standard deviation from fitting to 
R s (lo, T) spectra as a function of u> at each given T. The 
phenomenological parameter A _1 (T) is a measure of the 
energy-averaged quasiparticle transport relaxation rate, 
even though the energetics of the rf-wave superconducting 
state enter in a nontrivial manner ^J^MJ^ Ultimately, a 
full microscopic treatment of the charge transport should 
be used to fit the details of <7i(w,T). In the meantime, 
the phenomenological parameters provide useful insight. 

Upon examining the data for the Ortho-I sample in 
Fig-EH four key conclusions can be drawn. First, 1/A(T) 
appears to be a:b isotropic to within our experimental 
uncertainties. Second, fitting to the parameters arising 
from the resonator data indicates an average for the low 
temperature limit of A -1 ~ (45.8 ± 1.0) x 10 9 s" 1 and 
a temperature dependent component that grows rapidly, 
either as T x with x = 4.6±0.1, as indicated by the 
solid lines in the plot, or even exponentially. These ob- 
servations are consistent with those made in our previous 
publicationsii£*22i2£ Third, we note that the absolute val- 
ues of the damping 1/A(T) reported here (with y = 1.75) 
are essentially identical to the values of the Drude scat- 
tering rate t (T) (with y = 2) reported in our previ- 
ous work. We can therefore regard the parameter A in 
our phenomenological expression for the conductivity as 
a meaningful measure of the average quasiparticle trans- 
port lifetime. The fourth and final point is that sam- 
ple #1, with more complete oxygen filling, reveals lower 
1/A values, while the increased number of oxygen vacan- 
cies in sample #2 seem to cause a small broadening of 
the ai(u>) spectra. However, although the spectra are 
broader, the approximate isotropy of A is retained. 

Turning now to the Ortho-II results presented in Fig. [51 
a quite different scenario is observed. First, the spectral 
width parameter is evidently a:b anisotropic. Further- 
more, 1/A(T) remains strongly temperature-dependent 
all the way down to 1.2 K. This latter observation was 
a critical piece of evidence that led to the conclusion 
that the d-axis quasiparticle conductivity of Ortho-II or- 
dered YBa2Cu30g.5o is dominated by relatively weak (or 
Born-limit) scattering from static crystalline defects at 
low temperature^^ On the other hand, the linear-in- 
temperature regime apparently extrapolates back to a 
nonzero value of 1/A(T = 0), which is not in agreement 
with simple pictures of elastic scattering from weak point- 
like Coulomb potentials. Rather than scaling with T, 
we shoed that the width of the quasiparticle conductivi- 
tyspectrum scales as T + T with T = 2.0 KM 



The contrast between the overdoped and underdoped 
values of 1/A(T) is striking. As stated above, the spec- 
tral width gives the thermally-averaged quasiparticle re- 
laxation rate. Given that the chemical purity and crys- 
tallinity of all of the samples studied here are nominally 
the same, it is somewhat surprising that the conductiv- 
ity spectra for the two different dopings evolve in a very 
different manner. In fact, this observation suggests that 
oxygen doping in the CuO chain layer of YBa2Cu30g+ y 
has significant influence on the scattering of d-wave quasi- 
particles in the adjacent Cu02 planes. In addition, the 
large a:b anisotropy of the relaxation rate in Ortho-II 
is a surprising contrast to the approximate istotropy of 
Ortho-I. The anisotropy is puzzling since the low energy 
excitations reside near the four d-wave nodes along the di- 
agonals of the Brillouin zone, and these excitations must 
carry both the a and 6-axis currents. 



B. Spectral Weight 

Our phenomenological fits to the conductivity spectra 
yield the parameters <jo(T), A(T) and y(T). Rather than 
plot <7o(T), the zero- frequency intercept of the conduc- 
tivity, it is instead more insightful to present the inte- 
grated spectral weight of the quasiparticle conductivity. 
This quantity is often referred to as the normal fluid os- 
cillator strength. We obtain the spectral weight by inte- 
grating the fitted phenomenological form associated with 
the two-dimensional bands, and denote this quantity by 
(n n (T)e 2 /m*) 2D . 

The phenomenological spectrum in Eq. 1101 can be in- 
tegrated in closed form to give: 



n„(T)e s 



nr 



2D 



^ J dfla 2D (fl,T) 

2 v (T) 
y(T)A(T)sin(7r/y(T)) 



(14) 



Values for y(T) are given in Fig. H for the fits to bolo- 
metric measurements and values of y = 1.75 (sample #1, 
Ortho-I) and y = 1.45 (Ortho-II) are used for the fits 
to the cavity perturbation measurements. We use the 
values of A(T) plotted in Fig. By plotting the spec- 
tral weight we carry out an important check on both 
the experiments and the subsequent analysis of the data. 
The Ferrel-Tinkham-Glover sum rule^ states that the 
quasiparticle spectral weight and the superfluid spectral 
weight must add to a temperature-independent constant 
ne 2 /m*: 



E 



n n (T)e 2 



n s (T)e 2 



El ne 
\ TO* 



(15) 

Here, the London penetration depth of Eq.|8]is related to 
the superfluid spectral weight n s (T)e 2 /to* defined as the 
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sum over all bands that participate in superconductivity: 




It should be possible to observe the transfer of spec- 
tral weight between the quasiparticles and the supercon- 
ducting condensate as a function of temperature. The 
sum rule does not require complete condensation — it 
is entirely possible that residual quasiparticle oscillator 
strength will remain as T^O. We will demonstrate that 
in our YBa2Cu306+ y samples, up to 3% of the low fre- 
quency spectral weight fails to condense. A similar con- 
clusion has been drawn from microwave measurements 
on Tl2Ba2Cu06±<5 and Bi2Sr2CaCu20s, where at least 
10% of the oscillator strength remains as T— >0i22i 

In Figs. EI| an( l in) we show the integrated quasipar- 
ticle spectral weight from the two-dimensional bands of 
YBa2Cu306+y Representative error bars are shown at 
selected temperatures and were calculated from the stan- 
dard deviation of the best fit parameters. In addition, we 
show the corresponding loss of superfluid spectral weight, 
defined to be: 

n s (T = 0)e 2 n s (T)e 2 _ 1 1_ 

MoA 2 L (T = 0) moA|(T)' 1 ) 

for both dopings and for a and b directions. 

Errors in the choice of A^(T— >0) will have a strong 
effect on the absolute values of spectral weight obtained 
from our analysis. However, by fortunate coincidence, 
the uncertainty in Xl (T—*0) affects the the quasiparticle 
oscillator strength and the superfluid spectral weight in 
the same way. This can be seen by first noting that 
at low temperature and low frequency o\ oc R S /X L rs 
R S /X\{T=Q). Second, at low temperature one can write 
the following approximation for Eq. 1171 

n s (T = 0)e 2 n s (T)e 2 _ 2AA L (T) 

m* m* ~ // o Al(T=0)' 1 ' 

where AX L (T) = X L (T) - X L (T = 0). Therefore the 
change with temperature of the integrated quasiparticle 
spectral weight and the superfluid spectral weight both 
scale as 1/X\(T=Q). This is an important result because 
it allows us to test our data and analysis for consistency 
with the oscillator strength sum rule in a manner that is 
insensitive to uncertainties in Xl(T = 0). 

There are three features common to all panels in 
Figs. ^3 and El First, the superfluid spectral weight 
determined from measurements of AA(T) varies linearly 
with temperature at low T. This is a general conse- 
quence of having line nodes in the superconducting order 
parameter £ Second, the quasiparticle spectral weight at- 
tributed to the two-dimensional bands also varies linearly 
with temperature at low T. The low T slopes were ob- 
tained by linear fitting to the data of Figs.lllland ll2l and 
the results are tabulated in Table[H] Third, all of the in- 
tegrated quasiparticle spectral weight data extrapolate to 
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FIG. 11: Normal fluid oscillator strength for fully doped 
YBa2Cu30e.993 associated with the two-dimensional CUO2 
planar bands. Measurements of the real and imaginary parts 
of the conductivity allow this to be determined by two meth- 
ods: from integrating the 0-1(01) data (see Eg. 1101 and Fig.|^J, 
and from the disappearance of oscillator strength in the su- 
perfluid response, as determined by measurements of AA(T) 
at 1.1 GHz (see Eg. 1171 . The close agreement in slopes for 
the a-axis indicates that the conductivity sum rule is obeyed. 
This is not the case for the 6-axis where a considerable frac- 
tion of the total spectral weight resides in the broad spectral 
feature in eri(w) associated with the one-dimensional CuO 
chains. The difference in the fo-axis case between sample #1 
and sample #2 indicates that this parameter is very sensitive 
to oxygen content. 



nonzero values at T = indicating that a finite fraction of 
the spectral weight does not condense. Our estimates of 
the T — intercepts are tabulated in Table 11111 Approxi- 
mately 0.5% and 3% of the total spectral weight remains 
uncondensed for the Ortho-I YBa2Cu30g.993 and Ortho- 
II YBa2Cu3C>6.5o material respectively. The large error 
bars for the 6-axis results in Tables ITfl and Hill are due to 
the uncertainty in the fit parameter cr\ (u> — > 0,T) used 
to capture contributions from the CuO chain electrons. 

The results in Table [H] lead to a critical conclusion: 
along the a-axis, for both dopings, the rate at which the 
superfluid spectral weight decreases with temperature is 
equal in magnitude to the rate at which quasiparticle 
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FIG. 12: Normal fluid oscillator strength for Ortho-II ordered 
YBa2Cu3O6.50 associated with the two-dimensional CuCh 
planar bands. As for the fully doped case, the a-axis data 
closely obey the sum rule and the fc-axis data are more com- 
plicated. Note that only integration of o\ (uj) gives the correct 
T — > offset since our calculation of n„e 2 /m* from 1/A 2 as- 
sumes there is no residual quasiparticle oscillator strength. In 
all cases we observe a fraction of the total spectral weight that 
remains uncondensed (see Table ITTTt similar to observations 
in other cuprates. 



spectral weight increases with temperature. Our results 
are therefore consistent with the Ferrel-Tinkham-Glover 
sum rule in both Ortho-I YBa2Cu3 06.993 and Ortho-II 
YBa 2 Cu306.5o- This gives considerable support for the 
general correctness of the phenomenological lineshape of 
Eq. I1UI Again, we emphasize that this agreement is in- 
dependent of the exact choice of A£(T — > 0). 

The 6-axis results in Table [n] reveal that the rate at 
which the quasiparticle spectral weight attributed to the 
two-dimensional bands increases with temperature is less 
than the rate at which superfluid spectral weight de- 
creases with temperature for both dopings. This obser- 
vation is consistent with the idea that part of the quasi- 
particle spectral weight resides in the broad spectral fea- 
ture associated with the quasi-one-dimensional chain-like 
band. Previously^ we have shown that this is a very 
broad feature and therefore we cannot estimate its spec- 



Doping 


Orientation 


Spectral Weight 


Slope 








(lO^^msK)- 1 ) 


Ortho-I 


a 


2D-QP 


5.4 ±0.7 


sample #1 


a 


from 1/A 2 (T) 


5.4 




b 


2D-QP 


5.2 ±0.7 




b 


from 1/A 2 (T) 


7.9 


Ortho-I 


a 


2D-QP 


4.8 ±0.2 


sample #2 


b 


2D-QP 


5.9 ±0.1 


Ortho-II 


a 


2D-QP 


2.0 ±0.5 




a 


from 1/A 2 (T) 


2.0 




b 


2D-QP 


4.0 ± 1.8 




b 


from 1/A 2 (T) 


4.8 



TABLE II: Temperature slopes of the spectral weight 
n„(T)e 2 /m* from the data shown in Figs. llOl and llll The two- 
dimensional quasiparticle (2D-QP) slopes given in this table 
are an average of the slopes of the low temperature resonator 
and bolometry data in the cases where both measurements 
exist. 



tral weight. The fit parameter <j\ d {u> — > 0, T) for both 
dopings is plotted in Fig. 1131 Error bars are shown at 
selected temperatures and indicate one standard devia- 
tion in this fit parameter. The resonator measurements 
on Ortho-I suggest that a\ D {u — > 0,T) varies roughly 
linearly with T and extrapolates to a T = intercept of 
(7± 3) x 10 6 n _1 m _1 . The bolometry results on sam- 
ple #2 suggest a much lower T = intercept, although 
the error estimates are substantial. The differences be- 
tween sample #1 and sample #2 are most likely due to 
the different filling of the CuO chains in the two sam- 
ples. Sample #1 with T c = 89 K has an oxygen content 
of y = 0.993, while sample #2 with T c — 91 K was not 
so carefully controlled. The extent to which the off-plane 
oxygen atoms affect the microwave conductivity is un- 
clear and will be the subject of future studies. 

Extracting precise values of a\ D '(oj — > 0, T) from 
the Ortho-I bolometry data is difficult since the spec- 
tral width of the narrow component is a sizeable frac- 
tion of the experiment's bandwidth at all temperatures. 
In the Ortho-II case we observe much narrower low T 
spectra which allow more unambiguous extractions of 
o\ d {uj -> 0,T). For Ortho-II, a\ D {u -> 0,T) varies non- 
linearly with T and approaches a low T limiting value 
of (2± 1) x lO 6 ^ 1 !!!- 1 . Despite these differences, a\ D 
does not vanish as T — *■ for either Ortho-I or Ortho-II 
YP^CusOe+j/. This observation is consistent with the 
theoretical work of Atkinson— that suggests oxygen va- 
cancies in the CuO chains of YBa2Cu306+ 3/ may serve 
as strong pair breakers. 



VI. CONCLUSIONS 

The similarities and differences between the data 
shown for YBa 2 Cu306.5o and YBa 2 Cu3 06.993 offer in- 
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FIG. 13: Temperature dependence of the broad spectral fea- 
ture in the fe-axis conductivity associated with the quasi-one- 
dimensional CuO chain bands for both Ortho-I and Ortho-II 
YBa2Cu3C>6+i/. The conductivity component a\ D {uj — > 0, T) 
in Eq. 1131 is frequency independent in our phenomenological 
model since it is known to possess a width much greater than 
our highest measurement frequency. 



Doping 


Orientation 


RQOS 






(10 17 (fims)" 1 ) 


Ortho-I 


a 


6±2 


sample #1 


b 


10 ±7 


Ortho-I 


a 


7± 1 


sample #2 


b 


7.2 ± 0.6 


Ortho-II 


a 


6±2 




b 


10 ±6 



TABLE III: Residual quasiparticle oscillator strength 
(RQOS) from the bolometry data shown in Figs. HUl and 1111 
At most this corresponds to 3% of the total spectral weight. 



sight into the nature of the in-plane anisotropy in this 
system, the doping dependence of the inelastic scatter- 
ing, and the elastic scattering of nodal quasiparticles by 
defects. Both materials show a very broad component in 
<j{uj) in the direction parallel to the CuO chains, which 
is presumably due to the quasi-one-dimensional Fermi 
surface sheet derived from a CuO chain band hybridized 
with the Cu02 plane bands. Atkinson showed that such 
a mechanism is needed to explain the anisotropy of the 
in-plane superfluid density^ He also needed to invoke 
strong pair-breaking by O defects to explain the lack of 
a low temperature upturn in the 6-axis superfluid den- 
sity, an upturn expected if the superconductivity on the 



chains arises from the proximity effect. This strong scat- 
tering would likely explain the broad <t(uj) spectrum of 
quasiparticle excitations on this Fermi surface sheet. 

An important and puzzling issue arises by examining 
the anisotropy of the transport relaxation rate and tem- 
perature dependent spectral weight. In the cuprates, the 
Cu02 plane currents are carried by nodal quasiparticles. 
Since the gap nodes lie close to the diagonals of the Bril- 
loiun zone, these excitations carry both a and 6-axis cur- 
rents. We might therefore expect the transport relax- 
ation rate, as measured by the width of the <j\ d {uj,T) 
spectra, to be a—b plane isotropic. This is indeed the case 
in the Ortho-I sample: the conductivity spectra have ap- 
proximately the same width in the a and b directions, 
and the temperature-dependent spectral weights agree 
to within experimental uncertainty. The latter conclu- 
sion is sensitive to the choice of Xl{T — ► 0), but we 
conclude that Ortho-I YBa2Cu3 06.993 is well described 
by a 3-band model consisting of a CuO chain-like band 
acting in parallel with the bonding/antibonding Cu02 
plane bands. However, the Ortho-II results present a 
surprising contrast. The widths of the conductivity spec- 
tra and the temperature-dependent quasiparticle spec- 
tral weights differ by a factor of two between the a and 
b directions. One of the possible explanations for this 
anisotropy may be the effect of Ortho-II ordering on the 
electronic structure of the material. Ortho-II ordering 
doubles the unit cell of YBa2Cu30g+ y in the a-axis direc- 
tion, and can therefore substantially alter the Fermi sur- 
face topology with respect to the fully-doped compound. 
In particular, Bascones et al. have predicted that Ortho- 
II ordering can generate two new quasi-one-dimensional 
bands that are derived entirely from the Cu02 plane 
stated This might generate additional nodal quasiparti- 
cles derived from the new quasi-one-dimensional bands, 
increasing the a : b anisotropy of the in-plane conductiv- 
ity. We leave this issue as an open question for future 
work in this area, and as a feature that must be kept in 
mind for any theoretical treatment of this phase. 

Figure demonstrates that the rapid collapse of the 
inelastic scattering below T c occurs in both of these 
highly-ordered phases of YBa 2 Cu306+y One theoret- 
ical proposal treats the opening of a gap in the spec- 
trum of excitations that scatters the nodal quasiparti- 
cles, including the key detail that Umklapp processes are 
needed to relax charge currents . 18 i 19 The T dependence 
of the damping 1/A is consistent with exponentially- 
activated behaviour, the expectation for Umklapp scat- 
tering. The activated behaviour arises because there is a 
minimum threshold quasiparticle wave vector that satis- 
fies the Umklapp requirement. The new results for the 
Ortho-II phase YBa2Cu306.so are qualitatively consis- 
tent with this picture. The damping rises even more 
quickly, broadening ai(u>) out of the microwave range by 
20 K. This could result not only from the lower T c , but 
also because changes in the Fermi surface at this lower 
doping would enhance Umklapp processes by moving the 
required excitation closer to the nodes as compared to 
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YBa2Cu3 6.993- Unfortunately, a quantitative fit is dif- 
ficult because of the strong linear term in the elastic scat- 
tering in this sample. 

Quasiparticle scattering in YBa2Cu3 6.993 and 
YBa2Cu30e.5o samples is likely dominated by elastic 
scattering below about 20 and 10 K, respectively. 
In both cases, the lowest temperature spectrum at 
1.3 K has the cusp-like shape expected for weak 
scattering. 22 ! 50 ! 51 However, the YBa2Cu3C>6.5o spectra 
retain the weak-scattering cusp-like form up to higher 
temperatures where inelastic scattering takes over, but 
the YBa2Cu3 06.993 spectra take on a more Drudc like 
shape and do not show the temperature independent 
<7i (u> — > 0) expected for weak limit scattering. One 
possible explanation is that the scatterers are neither in 
the Born nor the unitary limit, something also suggested 
by recent thermal conductivity measurements^ In such 
a case, the shape and temperature evolution depend on 
the strength of the scattering, the density of scatterers 
and the size of the superconducting energy gap. Detailed 
fits to the spectra will be the subject of future work, 
but one can qualitatively state that the scattering in 
YBa2Cu3C>6.993 appears to be slightly farther from the 
Born limit than it is in YBa2Cu306.5o- It has also been 
suggested that the nearly linear temperature dependence 
of ai(u>,T) for any given fixed frequency u> > at low 
T is also a consequence of intermediate scattering phase 
shiftsi&S 

In attempting to generate a microscopic description 
of conductivity data many researchers have focused on 
models of quasiparticle scattering from point-like defects. 
This is a logical starting point, but it has become ap- 
parent that such models are incapable of simultaneously 
accounting for the observed spectral lineshapes and the 
significant amount of oscillator strength that remains un- 
condensed at 1.2 K, more than would be expected as the 
conductivity evolves towards the T = limit. Schachingcr 
et al. tried to invoke intermediate scattering phase shifts 
for point scatterers to explain experimental results, but 
did not present a quantitative set of parameters that 
could simultaneously describe the spectral shape, its 
temperature dependence, and the uncondensed oscilla- 
tor strength at 1.2 This problem is widespread and 



far more severe in other cuprates. In Bi2Sr2CaCu20s+5 
for instance there is much greater uncondensed oscilla- 
tor strength observed in both terahertz measurements of 
films^ and microwave measurements of single crystals^ 
In these cases, residual oscillator strength as high as 30 % 
of the superfluid density was observed, compared to the 
few percent seen in the YBa2Cu306+y samples discussed 
here. Motivated by the spectroscopic inhomogeneity ob- 
served in scanning tunneling spectroscopy^ OrensteinSi 
has suggested that this oscillator strength was due to a 
collective mode in the presence of inhomogeneity. More 
recently, Nunner et al. have been driven by tunneling 
results to suggest that randomly distributed interstitial 
O dopants may be the source of the electronic inhomo- 
geneity in Bi 2 Sr 2 CaCu208+5i 5 ^ This has suggested that 
rather than focussing on point scatterers, a better ap- 
proach is to study extended scatterers, which give rise 
to slowly varying potentials on the CuC>2 planes created 
by off-plane disorder. Such models start from knowledge 
of the actual defects in the system and have had con- 
siderable success in fitting both the Bi 2 Sr2CaCu208+5 
microwave data and the YBa2Cu3 06.993 data presented 
here. The basic suggestion is that off-plane disorder 
is a controlling feature of all cuprates and that even 
though the relative amount disorder is much smaller in 
the oxygen-ordered phases of YBa2Cu306+ y , it still plays 
an important role in understanding its electronic prop- 
erties. The new spectroscopic data presented here, to- 
gether with that on other materials 5 ^ 5 ^ provides a test- 
ing ground for the different models: collective modes in 
an inhomogeneous superconductor, point scatterers with 
intermediate phase shifts, and scattering by off-plane dis- 
order. 
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